The pressure of QED from the 
two-loop 2PI effective action 



Sz. Borsanyi* a and U. Reinosa^ 6 , 

a Department of Physics and Astronomy, University of Sussex, 
Brighton, East Sussex BN1 9QH, United Kingdom. 

b Institut fur Theoretische Physik, Universitat Heidelberg, 
Philosophenweg 16, 69120 Heidelberg, Germany. 



Abstract 

We compute the pressure of hot quantum electrodynamics from the 
two-loop truncation of the 2PI effective action. Since the 2PI resum- 
mation guarantees gauge-fixing independence only up to the order of 
the truncation, our result for the pressure presents a gauge dependent 
contribution of C(e 4 ). We numerically characterize the credibility of 
this gauge-dependent calculation and find that the uncertainty due to 
gauge parameter dependence is under control for £ < 1. Our calcu- 
lation also suggests that the choice of Landau gauge may minimize 
gauge-dependent effects. 
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The diagrammatic approach to finite temperature quantum field theories 
heavily relies on the convergence properties of the used expansion scheme. 
To cure the poor performance of the perturbative loop expansion, various 
resummation schemes have been invented Among these, the loop ex- 
pansion of the two-particle-irreducible (2PI) effective action implements a 
ladder resummation, which respects thermodynamical consistency and en- 
ergy conservation [21 E] . These features make the 2PI scheme attractive for 
nonequilibrium field theory applications [I] . A prerequisite for a nonequilib- 
rium method to be credible is, however, its reliability in equilibrium. There, 
it is important to check the convergence of expansion series of the 2PI effec- 
tive action. To this aim we calculated the notoriously ill-behaved pressure in 
a scalar context in Ref . [5] , and found a monotonous dependence on the cou- 
pling constant as well as a relatively small next-to-leading order correction 
to the pressure value even at couplings of 0(1). 

In the framework of gauge theories, however, the implementation of this 
approximation scheme suffers from various difficulties. One of these is that 
thermodynamic observables computed within this scheme are gauge fixing 
independent only up to the order of the truncation. One can illustrate this 
point by considering the issue of gauge parameter dependence in the co- 
variant gauge. For vanishing background fields, the 2PI effective action is 
a functional of the fermion, gauge and ghost propagators (respectively de- 
noted by D, G and G g h) which also depends on the gauge-fixing parameter £: 
F 2 pi[D, G, G g b, £]. The thermal pressure of the system is obtained by eval- 
uating r 2 pi at its stationary poin10 D = D, G — G, G g h = Ggh, for a given 
temperature T, and by subtracting the same calculation at zero temperature: 



V = -^T 2Pl [D,G,G gh ^] 



T=0 



It is possible to show that the ^-dependence of V uniquely comes from the 
explicit ^-dependence of r 2 pi and that it disappears if, in Fourier space, 

J2w»°M=t- ( 2 ) 

This last equation is the BRST identity for the propagator of the exact 
theory [6], which may break in a truncated resummation. Indeed, within a 
given truncation of the 2PI effective action, the gauge symmetry does not 
impose any constraint [7j on the two-point function G above the order of 



lr The barred propagators denote the solution of the stationarity equations: 
ST 2P i/SD = 0, ST 2P i/SG = and ST 2PI /SG gh = 0. 
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the truncation^ Thus beyond this order Eq. (^j breaks and ^-dependent 
contributions to the pressure appear^] 

A certain number of strategies can be put forward in order to try to 
cope with these inconvenient features. The first possibility is to introduce 
further approximations, on top of the loop expansion. This is the case of 
the approximately self- consistent resummations introduced in Ref. [S]. Using 
this method, a gauge independent determination of the entropy of QCD 
has been possible and shows a good agreement with lattice results down to 
temperatures about 2.5 times the transition temperature. There is however 
no general understanding on how to systematize this approach and evaluate 
higher orders in a gauge independent manner. 

Another possibility is to stick to the loop expansion of the 2PI effective 
action but play with the freedom in the choice of field representations. In- 
deed, the exact theory is invariant under reparametrization of the fields and 
one could exploit this feature in order to define a loop expansion obeying 
certain properties. This idea has been discussed in Ref. [10] where it has 
been applied to the linear sigma model in order to define a systematic loop 
expansion of the 2PI effective action fulfilling Goldstone's theorem at any 
order of approximation. Unfortunately no field representation is yet known 
in gauge theories which would ensure that the BRST identity (j2J) is fulfilled. 

It is finally possible to isolate gauge independent terms in the expression 
of the pressure by separating contributions from different perturbative or- 
ders. This means a re-expansion of the propagators D and G in powers of 
the coupling. The resulting modified resummation scheme did not show a 
substantial improvement of convergence [TT] . 

A different point of view is based on the experience that the 2PI loop 
expansion is known to have good convergence properties One can 
expect that contributions above the order of accuracy, and in particular gauge 
dependences are under control, at least in a large range of coupling values. In 
this paper we explore this possibility and consider the two-loop truncation of 
the 2PI effective action using the standard parametrization of the fields. We 
work in the covariant gauge with arbitrary gauge-fixing parameter £, which 
allows us to study how large gauge dependent contributions can be. 

Before embarking on a numerical evaluation, one has however to pay 
special attention to a second aspect, namely that of renormalization. The 
difficulty is related to the fact that truncations of the 2PI effective action 
only resum particular subclasses of perturbative diagrams for which (per- 

2 Below this order, the truncated two-point function G coincides with the exact one 
which fulfills the BRST identity. 

3 More precisely, if one truncates the 2PI effective action at L-loop order, one expects 
gauge dependences to appear at order e 2L [5]. 
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turbative) theorems do not apply. Recently a large effort has been put into 
extending renormalization theorems to the particular classes of diagrams re- 
summed by the loop expansion of the 2PI effective action. This has been 
first achieved in the framework of scalar theories [13] as well as scalar the- 
ories coupled to a fermionic field [Hj, and more recently in the framework 
of QED [T5] in the covariant gauge. In this latter case, it is important to 
emphasize that the renormalization procedure differs substantially from the 
one in perturbation theory. The reason for this is that, for a given loop 
truncation of the 2PI effective action and in contrast to what happens in 
perturbation theory, the photon two- and four-point functions develop lon- 
gitudinal quantum and thermal corrections. Although these contributions 
are formally of higher order than the order of the truncation, they bring UV 
divergences which need to be removed before defining a continuum limit. In 
Ref. pj?] a renormalization procedure involving a new class of counterterms 
has been put forward which allows one to deal with this new kind of UV 
divergences and thus opens the way to practical calculations. 

In this letter, we apply these ideas in order to evaluate the pressure of 
QED from the two-loop 2PI effective action. By solving the stationarity 
equations for the propagators G and D and properly determining the coun- 
terterms at the temperature of interest we calculate the pressure as given in 
Eq. [U This numerical procedure is repeated for several gauge fixing parame- 
ters and for various renormalization scales in order to explore how severe the 
problem of gauge dependence is in the context of thermodynamic calculations 
within the 2PI framework. 

In this work, where we discuss gauge parameter dependence, it is essential 
that the considered discretization respects gauge symmetry. In this way, the 
only source for gauge dependences is the particular truncation we use. For 
numerical purposes it is also convenient to use lattice rather than dimensional 
regularization. We thus consider QED on a hypercubic lattice of spacing a. 
We denote by Np the number of points on the time direction and iV the 
number of points on each of the spatial directions. The inverse temperature 
is P — Npa and the spatial volume V = N 3 a 3 . We decompose the lattice 
action in three pieces: S = S g + S g f + Sf. As gauge-field action, we consider 
the non-compact action 



where the field-strength tensor F^x) = Aj l A u (x) — AlA^(x) is expressed in 
terms of the forward derivativ^ A^A u (x) = a" 1 [A v {x + p,) — A v {x)\. We use 

4 The notation /t stands for the vector of length a along the positive ijl direction. 




(3) 
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a discretized covariant gauge fixing term 

s *f = 4 fl4 £ E A ^(*) A ^(*) . (4) 
a? /if 

given in terms of the backward derivative A^A„(a;) = a~ l [A ll (x) — A u (x — fi)] 
for latter convenience. Finally, the fermionic action is taken to be the naive 
chiral action 

Sf = ~y a 4 \^( x + P-hi^M^) ~ ^0)7/^0)^0 + p)] , (5) 

X 

where U^x) = exp (iaeA fl (x)) represents a link variable. 

Normally, the interacting two-point function D(x,y) or G(x,y) corre- 
sponds to the correlator of two operators at x and y. On the lattice however, 
where the fundamental objects are link variables, it is more convenient to 
introduce these two-point functions as 

D(x, y) = (ip(x)-ip(y)) c and G^(x, y) = (A^(x)A u (y - (1 - 5^)z>)) c . (6) 

This definition maintains the usual translation and reflection symmetries of 
G and D, as well as the identity G^ u (x,y) = G u ^(x,y). Notice also that the 
discretization we consider here, respects the chiral symmetry of our massless 
fermion: D(x,y) = ^ 7 Ai Z) M (x, y). We shall thus consider the 2PI effective 
action as a functional of rather than a functional of D. 

Since the pressure ([T]) cannot be determined exactly, we consider the loop 
expansion of the 2PI effective action as obtained from the Cornwall- Jackiw- 
Tomboulis formula [3J (a trivial term stemming from the ghosts is included 
in our numerics but not written explicitly here): 

T 2P1 [D,G] = -NjTt [\ogD- 1 +Do 1 D] + ± Tr [log Gq X G\ + r int [-D, G] 

(7) 

where we have defined TrO = a 4 J2 X Si Ou(x, x) = PV J2i^a( x = 0) an d 
we have included the possibility of an arbitrary number of fermionic flavors 
Nj^ The functional T mt [D, G] is given - up to an overall sign - by all 0-leg 
2PI diagrams that one can draw using the two-point functions D and G and 
the tree level vertices generated by the lattice action. These arise from the 

5 The parameter Nf will be used in what follows in order to eliminate the doublers 
which appear as a result of discretizing the fermionic action. Since our discretization 
generates 16 fermion tastes (one pair in each direction) we shall set Nf to 1/16. This is 
similar, for instance, to the fourth root taken on the staggered fermion determinant in the 
context of lattice gauge field theory [IB] . 
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expansion of the link variable U^x), and in turn of Sf, in powers of A^x). 
To make sure that the pressure we calculate is correct up to 0(e 3 ), we have to 
expand U^x) to 0(e 2 ). The vertex obtained from expanding U^(x) to C(e 3 ) 
brings no contribution to the pressure in the case of vanishing background 
fields, which we assume throughout this work. 

Combining the 0(e) and 0(e 2 ) vertices into two-loop 2PI diagrams, per- 
forming the relevant traces and making use of the properties of D, we obtain 
the following contributions to the interacting part r; nt of the 2PI effective 
action: 

-^r mt = e 2 N f a 4 ^G llu {x)\D ll {x)D v {x + fi + P)+D v {x)D tM {x + p + i>) 

X, [i^v 

+ D^(x + u)D u (x + p) + D v (x + 0)D ll (x + p) 
+ e 2 N f a 4 £ G^(x) [2D^x - p)D^x + p) + 2D^x)D^x) 
- [Dy{x - p)D v {x + p) + D v {x)D v {x)\ 

V 

^r mt = a e 2 N f G m {x = 0) [d,(x = p) - D,(x = -p)] . (9) 

The contribution T mt is the usual fermion loop with a somewhat peculiar 
photon line (El). The lattice spacing a in T mt manifests that this diagram is a 
lattice artefact. Both fermion loops are individually quadratically divergent, 
they together make sure that at the lowest perturbative level the photon 
receives no mass renormalization. 

Together with a counterterm contribution 5Ti nt (see below), the expres- 
sions (EJ) and (IHl) provide the full 0(e 2 ) interaction part of the 2PI effective 
action: r- mt = T mt + T mt + 5Ti nt . If we now introduce the self-energies 

E,(x) = D-\x) - D'Kx) and U^(x) = G^(x) - G^{x) , (10) 

and use the explicit formula (|7j) for the 2PI effective action, we can write the 
stationarity equations defining the interacting two-point functions D and G 

as 



AN f E^x) = } and 11^ (x) = ^ J^T s ■ { 11 : 

a A pV oD^(x) a A pV oG^ u (x) 



The interacting two-point functions D and G are thus obtained after simul- 
taneously solving Eqs. (TlOl - flTT]) . As it can easily be checked, in the two-loop 
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approximation that we consider here, Eqs. do not involve any discretized 
integral in direct space. On the other hand, Eqs. ffTUj) can be conveniently 
solved in momentum space. To this order of the truncation we can thus 
completely avoid calculating loops by simply Fourier transforming the prop- 
agators back and forth in every step of the iterative procedure. We define 
the Fourier transforms of a generic fermionic (D) or gauge (G) two-point 
function, respectively, as 

r x D^k) = a 4 J2 e ' ik ' X D ^ x ) and ^(QG^k) =a 4 J2 G ^ k ) > 

X X 

(12) 

where ot^ik) — 1 if \i = v and a^ u (k) = exp(—ia(k^ + k u )/2) otherwise. This 
particular definition of the Fourier transform of G is connected to the fact 
that the gauge field has to be thought as attached to the midpoints of the 
links. Even this unusual variant of the fast Fourier transformation is available 
as legacy code [IT]. Solving Eqs. ffTUj) in Fourier space, needs that we deter- 
mine the Fourier transforms of the free inverse propagators. After inspection 
of the free (quadratic) contribution to S, one obtains 

Do,l(k) = ~h ^d G ^(k) = k 2 5^-(l-C 1 )kk»i (13) 

with the usual short-hand notations k^a = sm(k^a) and k^a = 2sin(k fl a/2). 

Eqs. ( II 01) -( fill) can be solved for any non- vanishing lattice spacing a, lead- 
ing to perfectly finite two-point functions D and G. In order to define a 
proper continuum limit of the latter, as a — > 0, one needs however to absorb 
UV divergences. Renormalization of Eqs. (1101) and (1111) was considered in 
Ref. [IS] in the context of dimensional regularization and at zero tempera- 
ture. There, renormalization was achieved by adding a contribution <5Fi nt to 
the functional T^t. This contribution carries the counterterms needed for 
renormalization. In extending this result to lattice regularization, one has to 
pay attention to the presence of new vertices originating from the expansion 
of the link variable U^(x) in powers of the field A^(x) (see above). In our 
present calculation, in addition to the usual vertex coupling A to ijj and ip 
(which leads to Tf nt ), there is a new vertex coupling A 2 to ip and ip (which 
leads to rf nt ). This new vertex brings an extra factor of a, which is such 
that the superficial degree of divergence of a given diagram is the same as in 
dimensional regularizationj^ It follows that we can here apply the same type 
of analysis of UV divergences as the one used in Ref. [T3] . 

6 Namely <5 = 4 — Ea — (3/2)£"^, where Ea and E^ respectively denote the number of 
external photon and fcrmion legs of the diagram at hand. 
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At two-loop order, the shift 5r int is given in lattice regularization by 



k,fj, 



q,v 



4 pV 



+ 



fiu k g 



+ 5M 2 5, 



k fj.v 

-4N f 5Z 2 J2h D »( k ) 

k,n 



(14) 



It leads to additional contributions at the level of the self-energies, in partic- 
ular a longitudinal wave function renormalization (8X) as well as a photon 
mass counterterm (6M 2 ). 

The counterterms 5gi and 5g 2 allow to remove sub divergences hidden in 
Eqs. (II II) and involving four photon legs (see below). After these have been 
removed, there only remain temperature independent overall divergences that 
need to be absorbed in the counterterms 6Z 2 , SZ 3 , 5X and 5M 2 . Although 
the exact 0(4) symmetry is broken on the lattice, the tensor structure of 
the self energies at a fixed scale k is restored in the continuum limit of an 
isotropic lattice theory. This allows us to use the renormalization conditions 
introduced in the context of the continuum theory [18]. In particular, one 
can show that the overall divergences have the structure 

^T(k) = -ak^ and Tl^(k) = n M 5, u + M^uk 2 - k,K) + ^k u (15) 

where a, n?, t^l an d t^m represent quantities which diverge as a — > 
(quadratically for ttm and logarithmically for the rest of them). Comparing 
these expressions to those for the counterterms, we find that all divergences 
can be absorbed by setting 



5Z 2 = —cr , 5Z- 3 = —txt , 5\ = —txl and 5M 2 = —tvm 



(16) 



The set of Eqs. (1161) does not fix the finite parts of the counterterms. In 
order to do so, we fix 5Z 2 , 5Z 3 , 5X and 5M 2 through the renormalization 
conditions 



3k* 



on 



22 



dk 2 



33 



dk 2 



, and II 



33 I A;* 



(17) 



where k* = (0, 0, fi, 0) and \x denotes our renormalization scale. The star on 
the self-energies means that these are considered at a reference temperature 
T*. The first two renormalization conditions are similar to those which are 
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used in perturbation theory and completely determine the counterterms 5Z 2 
and 5Z 3 . In perturbation theory, where the (lattice) Ward identity for TL(k) 



prevents the appearance of longitudinal corrections to the self energy, the 
third and fourth conditions in Eq. ( ITTj) are trivially satisfied. In our case, 
however, we need to fix two counterterms (5 A and 5M 2 ) that cancel UV 
divergences of 0(e A ). A natural way to fix these is to impose the Ward 
identity on IT at the renormalization point k* = (0,0, /i, 0) and for a given 
temperature T*. We do so at k* and in a small neighborhood of k*. In this 
way, we obtain the third and fourth renormalization conditions in Eq. ( TTTj) . 
The arbitrariness of this condition introduces an ambiguity of order 0(e 4 ). 

As already discussed in Ref . [15] , when renormalizing the two-point func- 
tion G, one has not only to pay attention to longitudinal overall diver- 
gences but also to longitudinal sub divergences which involve four-photon 
legs. Again, if no truncation is considered, these sub divergences automat- 
ically cancel since they reproduce the exact four-photon function which is 
transverse. However, for a given truncation of the 2PI effective action, this 
cancellation of divergences is only true up to the order of the truncation. 
Above, new divergences appear which need to be absorbed by means of the 
counterterms Sgi and Sg 2 - The particular structure of these divergences has 
been worked out in Ref. [15] for the case of dimensional regularization. The 
result is that, in order to absorb the four-photon divergences, one needs to 
impose, at the renormalization point, the transversality of a four-point func- 
tion defined by means of a set of Bethe-Salpeter equations. Here, we extend 
this result to the case of lattice regularization. 

The Bethe-Salpeter equations can be written as a closed set of equations 
for a four-point function V^ ap (p,k) involving four photon legs and a four- 
point function Wij tUp (p, k) involving two photon and two fermion legs [T5] . 
Similarly to what we did with the propagator D, we turn the Dirac indices i, j 
into one Lorentz index /x: W^apj^ = W^^p. Given that k* = (0, 0, fi, 0), 
the renormalization conditions fixing 5gi and 5g 2 , as given in Ref. [15], read 



In order to impose these renormalization conditions, we do not need to 
solve the set of Bethe-Salpeter equations for arbitrary values of the mo- 
menta and arbitrary configurations of Lorentz indices. Indeed, the set of 
equations remains closed if we fix one of the momenta to k = k* and two 




(18) 



VWfc*, k *) = and ^3333^, k*) = . 



(19) 
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of the Lorentz indices to a = p = 3. We shall thus consider equations for 
^V(p) = V^miVi k*) an d W a (p) = W at33 (p, k*). Introducing the notations 

A ap (p) = 6(p-k*)5 a3 5 p3 , (20) 
V IUf (p) = G lia (p)V a p(p)Gp v (p), (21) 

W» = -2 D» £ W,(p)5,(p) + W,(p) E D P (p)D p (p) , (22) 
p p 
we may write the corresponding set of Bethe-Salpeter equation as 

VM = --f 5 ^ E M?) ~ 2A fM)] - |j E - 2iV(g)] 



gMg) ^ ( a) (23) 
dD p {q) Wp{q) ' {26) 



5, per p v ' <j,p f \ -w 

We solved this pair of equations iteratively by adjusting Sgi and 8g2 after 
each step so that Eq. (fT9l) is always fulfilled. As expected, the numerical val- 
ues of these counterterms scale as ~ e 4 log(a). Once the counterterms have 
been fixed according to the renormalization conditions i fTTj) and (Tl9l) . we can 
solve for the physical two-point functions D and G which admit a proper 
continuum limit. Plugging this values into the CJT formula ([7]) truncated at 
two-loop order gives us a non-perturbative approximation to the QED pres- 
sure, compatible with perturbation theory up to order C(e 3 ). Notice that, 
even with all our counterterms, there is a quartic divergence remaining in the 
pressure. This divergence is temperature independent and can be removed 
by a 'cosmo logical constant' renormalization. The renormalization condition 
is usually given by the requirement of zero vacuum pressure. Here we do not 
renormalize or evaluate the model at zero temperature. We determine the 
counterterms in the equations of motion at T*. Then, using these countert- 
erms we evaluate the pressure at T* and T*/2. Assuming a ~ T 4 scaling 
with the temperature, we determine the pressure as the difference of the di- 
vergent pressure values as obtained from the formula of the effective action, 
divided by (15/16)(T*) 4 . The assumed scaling of temperature is broken due 
to the presence of the renormalization scale. This effect introduces an error 
of 0{e i ) which is above the actual accuracy of our calculation. 

In order to improve numerical stability, we took into account the fol- 
lowing points. Calculating the pressure difference involves the subtraction 
of two quartically divergent contributions. Instead, we carried out the spa- 
tial part of the trace in I^pi after performing the subtraction. An other 
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0.84 
0.82 

0.5 1 1.5 2 2.5 

e 

Figure 1: Two-loop QED pressure as a function of the coupling e and for 
different values of the gauge-fixing parameter £. The plain line corresponds 
to £ = (Landau gauge), long-dashed lines to £ = 1 (Feynman gauge) and 
short-dashed lines to £ = 2. The sensitivity with respect to the renormaliza- 
tion scale \i is illustrated in the case of the Feynman gauge. We also plot the 
perturbative C(e 2 ) result for comparison. 



important alteration to the equations above was the exclusion of the spa- 
tially homogeneous lattice mode on the level of the 2PI effective action. This 
was necessary to avoid instabilities as e — > 0, since the finite photon mass 
contribution behaves as ~ e 4 . 

In Fig. [H we plot the QED pressure in the two-loop 2PI approximation, 
for a wide range of coupling values (0 < e < 2.4) and for various values of the 
gauge-fixing parameter. As discussed in Ref. [H] the higher the gauge-fixing 
parameter is, the less convergent the 2PI loop expansion becomes. It is thus 
meaningless to consider our calculation for too high values of £ and, as sug- 
gested in Ref. [8], we restricted our calculations to values of the gauge- fixing 
parameter ranging from £ = (Landau gauge) to £ = 2. For small values of 
the coupling, our results are almost insensitive to the gauge fixing parame- 
ter and nicely reproduce the perturbative result to order 0(e 2 ). This comes 
as no surprise since the two-loop 2PI approximation contains all diagrams 
contributing to order (9(e 3 )Jll Numerically we find a good agreement with 
perturbation theory up to e ~ 1 which is precisely where the perturbative 

7 Our calculation reproduces the 0(e 3 ) result in the infinite volume limit only. 
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Figure 2: Renormalization scale dependence of the two-loop pressure. 



expansion usually breaks down. 

For large values of the coupling, our calculation becomes a priori sensitive 
to two types of uncertainties. First of all, renormalization is done by impos- 
ing renormalization conditions at a certain momentum k* = (0, 0, /i, 0) which 
introduces an artificial dependence on the scale /i. Moreover the truncation 
of the 2PI effective action introduces gauge parameter dependences starting 
at order 0(e A ). These two types of uncertainties can be taken as a way to 
estimate the error of the calculation. 

The dependence with respect to the renormalization scale /i is illustrated 
in Fig. [U for the case £ = 1 (Feynman gauge) where fi is varied in the interval 
7rT < fi < AttT as it is usually done in calculations at finite temperature. A 
study of the /i-dependence as the gauge-fixing parameter is varied and for a 
given value of the coupling is depicted in Fig. [2j Notice that, at fixed gauge- 
fixing parameter £, the /i-dependence is not monotonous. However /i = 2nT 
roughly represents the value at which the pressure reaches it maximum value, 
in this range. We notice that the uncertainty due to scale dependence is not 
particularly severe, which indicates the good convergence behavior of the 
2PI approach. Moreover this uncertainty is ~ 1% for £ = 2 and decreases 
considerably down to its minimum value reached for £ = 0, which makes the 
Landau gauge a particularly interesting choice among all possible gauges. We 
also notice that, in general, choosing a higher renormalization scale flattens 
the gauge dependence towards the Landau gauge value. 

The second source of uncertainties is gauge dependence. As already men- 
tioned a calculation for high values of the gauge fixing parameter makes little 
sense. In the considered range of gauge parameter values, the error due to 
gauge dependence is of the order of or less than 1 — 1.5%. The Landau gauge 
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plays again a special role since it corresponds to the value of £ for which 
the pressure is the less sensitive to gauge parameter dependence. Indeed, 
independently of the value of the coupling, one has — p^=o ~ £ 2 as £ — > 0, 
as it is clear on the logarithmic plot of Fig. [3j 




In conclusion, our calculation shows, in covariant gauge, a relatively small 
error coming from gauge parameter dependence. The parametric suppression 
of the gauge parameter dependence has already been shown in Ref. [5]. We 
have now established that the so far unknown coefficients of this parametric 
dependence do not spoil this behavior. The gauge dependence can also be 
regarded as a bonus feature, which opens a way to error estimates without 
the need for considering higher order diagrams. We think, that the 2PI 
effective action can be regarded as an efficient resummation technique for 
gauge theories, where the actual choice of gauge fixing has an impact on the 
quality of the resummation. As for the particular calculation presented here, 
the Landau gauge is the preferred choice. If this result persists in QCD, it 
could serve as a justification for the exclusive use of Landau gauge in the 
QCD Schwinger- Dyson equations [19J. 
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